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1 Newton's Laws: First principles 



The fundamental idea of kinematics is the discussion of the movement of objects, without actu- 
ally taking into account what caused the movement to occur. By using simple calculus, we can 
find all of the equations for kinematics. To simplify the learning process, we will only consider 
objects that move with constant acceleration. For the first few parts, we will also assume that 
there is no friction or air resistance acting on the objects. 
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2 Straight Line Motion (SLM) 



The name of this section Straight Line Motion means that we begin learning about the subject 
of kinematics by observing motion in one dimension. This means that we will only take one axis 
of a 3D (x, y, z) coordinate system into account. We will be using the x axis as our axis of motion. 

Throughout our discussion, we will look at the motion of a rigid body, one that does not deform 
as it moves. We idealize the rigid body by assuming that it has no dimensions and is infinitely 
small. This way we can talk about the entire body instead of saying, "The front of the body is at 
this point while the rear of the body is at some other point". Also, subscripts on the variables in 
our equations will indicate the initial value, i, and the final value, /, of that variable. 

2.1 Displacement 

To start off with we will define the term displacement. Displacement is basically the short route 
of getting from one place to another, a straight line from the starting point of motion to the 
ending point of motion. No matter how much movement takes place in between, we only care 
about the first location and the second location. We will use the variable, x, to stand for the 
locations of the rigid body that we are discussing. 

In order to define motion we first must be able to say how far an object has moved. This is 
done by subtracting the final value of the displacement by the initial value of the displacement. 
Or, in other words, we subtract the initial position of the object on our coordinate system from 
the final position of the object on our coordinate system. It is necessary to understand that the 
values which you receive for this variable can be either positive or negative depending on how 
the object is moving and where your coordinates start: 

X — ^ f "^i 

2.2 Velocity 

The term velocity, v, is often mistaken as being equivalent to the term speed. The term velocity 
refers to the displacement that an object traveled divided by the amount of time it took to move 
to its new coordinate. From the above discussion you can see that if the object moves backwards 
with respect to our coordinate system, then we will get a negative displacement. Since we cannot 
have a negative value for time, we will then get a negative value for velocity. The term speed, on 
the other hand refers to the magnitude of the velocity, so that it can only be a positive value. We 
will not be considering speed in this discussion, only velocity. By using our definition of velocity 
and the definition of displacement from above, we can express velocity mathematically like this: 
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Straight Line Motion (SLM) 



{Xf-Xi) 

v- — 

itf-ti) 

The line over the velocity means that you are finding the average velocity, not the velocity at a 
specific point. This equation can be rearranged in a variety of ways in order to solve problems 
in physics dealing with SLM. Also, by having acceleration as a constant value, we can find the 
average velocity of an object if we are given the initial and final values of an object's velocity: 

{Vf+Vi) 

v= — 

2 

These equations can be combined with the other equations to give useful relationships in order 
to solve straight line motion problems. For example, to find the initial position of an object, if we 
are given its final position and the times that the object began moving and finished moving, we 
can rearrange the equation like this: 

Xi = Xf - V(tf - ti) 

2.3 Acceleration 

The term acceleration means the change in velocity of an object per unit time. We use the vari- 
able, a, to stand for acceleration. Basically, the sign of acceleration tells us whether the velocity 
is increasing or decreasing, and its magnitude tells us how much the velocity is changing. In or- 
der for an object that is initially at rest to move, it needs to accelerate to a certain speed. During 
this acceleration, the object moves at a certain velocity at a specific time, and travels a certain 
distance within that time. Thus, we can describe acceleration mathematically like this: 

(Vf-Vi) 

a= — 

itf-ti) 

Again, we can rearrange our formula, this time using our definition for displacement and veloc- 
ity, to get a very useful relationship: 

1 2 
Xf = Xi + Viitf - ti) + -a{tf - t{) 

2.4 History of Dynamics 
2.4.1 Aristotle 

Aristotle expounded a view of dynamics which agrees closely with our everyday experience of 
the world. Objects only move when a force is exerted upon them. As soon as the force goes away, 
the object stops moving. The act of pushing a box across the floor illustrates this principle -- the 
box certainly doesn't move by itself! 
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Newton's 3 Laws of Motion 



However, if we try using Aristotle's dynamics to predict motion we soon run into problems. It 
suggests that objects under a constant force move with a fixed velocity but while gravity definitely 
feels like a constant force it clearly doesn't make objects move with constant velocity. A thrown 
ball can even reverse direction, under the influence of gravity alone. 

Eventually, people started looking for a view of dynamics that actually worked. Newton found 
the answer, partially inspired by the heavens. 

2.4.2 Newton 

In contrast to earthly behavior, the motions of celestial objects seem effortless. No obvious forces 
act to keep the planets in motion around the sun. In fact, it appears that celestial objects simply 
coast along at constant velocity unless something acts on them. 

This Newtonian view of dynamics — objects change their velocity rather than their position 
when a force is exerted on them — is expressed by Newton's second law: 

F = ma 

where F is the force exerted on a body, m is its mass, and a is its acceleration. Newton's first 
law, which states that an object remains at rest or in uniform motion unless a force acts on it, is 
actually a special case of Newton's second law which applies when F=0. 

It is no wonder that the first successes of Newtonian mechanics were in the celestial realm, 
namely in the predictions of planetary orbits. It took Newton's genius to realize that the same 
principles which guided the planets also applied to the earthly realm as well. 

In the Newtonian view, the tendency of objects to stop when we stop pushing on them is simply 
a consequence of frictional forces opposing the motion. Friction, which is so important on the 
earth, is negligible for planetary motions, which is why Newtonian dynamics is more obviously 
valid for celestial bodies. 

Note that the principle of relativity is closely related to Newtonian physics and is incompatible 
with pre-Newtonian views. After all, two reference frames moving relative to each other cannot 
be equivalent in the pre-Newtonian view, because objects with nothing pushing on them can 
only come to rest in one of the two reference frames! 

Einstein's relativity is often viewed as a repudiation of Newton, but this is far from the truth — 
Newtonian physics makes the theory of relativity possible through its invention of the principle 
of relativity. Compared with the differences between pre-Newtonian and Newtonian dynamics, 
the changes needed to go from Newtonian to Einsteinian physics constitute minor tinkering. 

2.5 Newton's 3 Laws of Motion 

2.6 Newton's first law: 

An object at rest tends to stay at rest and an object in motion tends to stay in motion with the 
same speed and in the same direction unless acted upon by an unbalanced force. 
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Straight Line Motion (SLM) 



2.7 Newton's second law: 

The acceleration of an object as produced by a net force is directly proportional to the magnitude 
of the net force, in the same direction as the net force, and inversely proportional to the mass of 
the object. 

2.8 Newton's third law: 

If body A exherts a force on body B, then body B exherts a force on body A that is equal in mag- 
nitude but opposite in direction from the force from A. 



When a force is exerted on an object, energy is transferred to the object. The amount of energy 
transferred is called the work done on the object. However, energy is only transferred if the ob- 
ject moves. Work can be thought of as the process of transforming energy from one form into 
another. The work W done is 



where the distance moved by the object is A x and the force exerted on it is F. Notice that work 
can either be positive or negative. The work is positive if the object being acted upon moves in 
the same direction as the force, with negative work occurring if the object moves opposite to the 
force. 

This equation assumes that the force remains constant over the full displacement. If it is not, 
then it is necessary to break up the displacement into a number of smaller displacements, over 
each of which the force can be assumed to be constant. The total work is then the sum of the 
works associated with each small displacement. In the infinitesimal limit this becomes an inte- 
gral 



If more than one force acts on an object, the works due to the different forces each add or subtract 
energy, depending on whether they are positive or negative. The total work is the sum of these 
individual works. 

There are two special cases in which the work done on an object is related to other quantities. 
If F is the total force acting on the object, then by Newton's second law W= FA x= mA x- a . 
However, o= dvl dt where v is the velocity of the object, and A x~vA t, where A t is the time 
required by the object to move through distance A x . The approximation becomes exact when 
A x and A t become very small. Putting all of this together results in 



2.9 Work 



W = FAx 
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Energy 



d mv 2 



W to t a i = m—vAt = At 



t 



dt 2 



We call the quantity mi//2 the kinetic energy, or K. It represents the amount of work stored as 
motion. We can then say 



Thus, when F is the only force, the total work on the object equals the change in kinetic energy 
of the object. 

The other special case occurs when the force depends only on position, but is not necessarily the 
total force acting on the object. In this case we can define a function 



and the work done by the force in moving from jq to %2 is V{ xi)- V{ X2), no matter how quickly 
or slowly the object moved. 

If the force is like this it is called conservative and V is called the potential energy. Differentiating 
the definition gives 




The minus sign in these equations is purely conventional. 
If a force is conservative, we can write the work done by it as 



where is the change in the potential energy of the object associated with the force of interest. 

2.10 Energy 

The sum of the potential and kinetic energies is constant. We call this constant the total energy 
E: 



If all the forces are involved are conservative we can equate this with the previous expression for 
work to get the following relationship between work, kinetic energy, and potential energy: 



W total = —At = AK 
at 




W= Ax= -AV 

dx 



E = K+U 



AK=W = -AU 
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Straight Line Motion (SLM) 



Following this, we have a very important formula, called the Conservation of Energy : 



MK+U) = 

This theorem states that the total amount of energy in a system is constant, and that energy can 
neither be created nor destroyed. 



2.11 Power 

The power associated with a force is simply the amount of work done by the force divided by 
the time interval over which it is done. It is therefore the energy per unit time transferred to the 
object by the force of interest. From above we see that the power is 

FAx 

P= = Fv 

At 

where is the velocity at which the object is moving. The total power is just the sum of the powers 
associated with each force. It equals the time rate of change of kinetic energy of the object: 

_ Wtotal _ dK 
Ftotal ~ At ~ dt 

We will now take a break from physics, and discuss the topics of partial derivatives. Further infor- 
mation about this topic can be found at the Partial Differential Section 1 in the Calculus 2 
book. 



2.12 Partial Derivatives 

In one dimension, the slope of a function, /( x), is described by a single number, df/dx. 

In higher dimensions, the slope depends on the direction. For example, if /= x+2 y, moving one 
unit x-ward increases / by 1 so the slope in the x direction is 1, but moving one unit y-ward 
increases / by 2 so the slope in the y direction is 2. 

It turns out that we can describe the slope in n dimensions with just n numbers, the partial 
derivatives of / 

To calculate them, we differentiate with respect to one coordinate, while holding all the others 
constant. They are written using a d rather than d. E.g. 

df 

f(x + dx,y,z) = f{x,y,z) + —dx (1) 



1 http://en.wikibooks.org/wiki/Calculus%2FPartial%20differential%20equations 

2 http://en.wikibooks.org/wiki/Calculus 
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Partial Derivatives 



Notice this is almost the same as the definition of the ordinary derivative. 

If we move a small distance in each direction, we can combine three equations like 1 to get 

df df df 

df = —dx+ — dy + — dz 
dx dy dz 

The change in / after a small displacement is the dot product of the displacement and a special 
vector 



f/^.{dx,dy,dz)= % ^f^f 



This vector is called the gradient of / It points up the direction of steepest slope. We will be 
using this vector quite frequentiy. 

2.12.1 Partial Derivatives #2 

Another way to approach differentiation of multiple -variable functions can be found in Feynman 
Lectures on Physics vol. 2. It's like this: 

The differentiation operator is defined like this: df = f{x + Ax, y + Ay, z + Az) - fix, y, z) in the 
limit of Ax, Ay, Az -» 0. Adding & subtracting some terms, we get 

df = (/(x+ Ax, y+ Ay, z+ Az) -fix, y+ Ay, z+ Az)) + {fix, y+ Ay, z+ Az) -fix, y, z+ Az)) + {fix, y, z+ 
Az)-f{x,y,z)) 



and this can also be written as 



df df df 
df = — dx + — dy + — dz 
dx dy dz 



2.12.2 Alternate Notations 



For simplicity, we will often use various standard abbreviations, so we can write most of the 
formulae on one line. This can make it easier to see the important details. 

We can abbreviate partial differentials with a subscript, e.g., 

D x h{x,y) = — D x D y h = D y D x h 

or 

dxh = — 
dx 

Mostly, to make the formulae even more compact, we will put the subscript on the function itself. 
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D x h = h x h xy = h yx 

2.13 More Info 

Refer to the Partial Differential Section 3 in the Calculus 4 book for more information. 

2.14 Motion in 2 and 3 Directions 

Previously, we discussed Newtonian dynamics in one dimension. Now that we are familiar with 
both vectors and partial differentiation, we can extend that discussion to two or three dimen- 
sions. 

Work becomes a dot product 

W = F-Ax 

and likewise power 

P = F-v 

If the force is at right angles to the direction of motion, no work will be done. 

In one dimension, we said a force was conservative if it was a function of position alone, or 
equivalently, the negative slope of a potential energy. 

The second definition extends to 

F= -W 

In two or more dimensions, these are not equivalent statements. To see this, consider 

Dyf X - D X Fy = V X y - Vy X 

Since it doesn't matter which order derivatives are taken in, the left hand side of this equation 
must be zero for any force which can be written as a gradient, but for an arbitrary force, depend- 
ing only on position, such as F=( y, - x, 0), the left hand side isn't zero. 

Conservative forces are useful because the total work done by them depends only on the differ- 
ence in potential energy at the endpoints, not on the path taken, from which the conservation of 
energy immediately follows. 

If this is the case, the work done by an infinitesimal displacement dx must be 



3 http://en.wikibooks.org/wiki/Calculus%2FPartial%20differential%20equations 

4 http://en.wikibooks.org/wiki/Calculus 
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Circular Motion 



W = V{x) - V(x + dx) = -(W) • dx 

Comparing this with the first equation above, we see that if we have a potential energy then we 
must have 



F = -W 

Any such F is a conservative force. 

2.15 Circular Motion 

An important example of motion in two dimensions is circular motion. 
Consider a mass, m, moving in a circle, radius r. 

The angular velocity, a) is the rate of change of angle with time. In time A t the mass moves 
through an angle A0= wA t. The distance the mass moves is then r sin A0, but this is approxi- 
mately r/S.6 for small angles. 

Thus, the distance moved in a small time A t is rwA t, and divided by A t gives us the speed, v. 



v = cor 

This is the speed not the velocity because it is not a vector. The velocity is a vector, with magni- 
tude a) r which points tangentially to the circle. 

The magnitude of the velocity is constant but its direction changes so the mass is being acceler- 
ated. 

By a similar argument to that above it can be shown that the magnitude of the acceleration is 



a = a)v 

and that it is pointed inwards, along the radius vector. This is called centripetal acceleration. 
By eliminating v or a) from these two equations we can write 

2 V\ 

a = -a) r = r 

r 

Year 10 Force and Work 7 periods 

• Possible Approach: 

• Lesson Suggested Learning Outcomes Course Outline Resources 

• 1-3 Recognise that unbalanced forces cause acceleration. 

• Investigate the relationship between the unbalanced force acting on an object, is mass *and 
acceleration: Newton's Second Law. (Use Fnet = ma giving correct units for all quantities) 
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•**Reinforce ideas about force causing acceleration (= the rate at which speed changes) 
• 

• Experiment 1: Finding the force required to move an object (static friction) 

• Discuss the factors that may affect the frictional forces acting on the block (e.g. load on *the 
block, surface area in contact with the bench, type of the surfaces in contact . . . ) 

• Attach a force meter to a block at rest on the bench. Pull the force meter and slowly increase 
*the force until the block just begins to move - record the force meter reading, tabulate data, 
"repeat, calculate average. 

• [Or, use a bag to which washers may be added, attached to a cord and hanging over the end of 
*the bench to pull the block.] 

• quantities and the third unknown) - a three symbol triangle may simplify the maths. 
•••••• Homework: Coursebook- Newton's second Law 

• force (e.g. when a force acts on a body, it moves faster and kinetic energy increases) 

• • Introduce the relationship between energy transferred and work done. 

• • Define work done in terms of the applied force and the distance through which 
this ********force acts (Work done = force x distance) 

....... u se the work rule to calculate work / force / distance (given two quantities and the 

*third unknown) - a three symbol triangle may simplify the maths 



Homework: Coursebook - Work Done and Energy 

• 6 - 7 Investigate and describe the pattern of results formed from graphing 
*the effects of applied forces on a spring. 



• Explore the relationships between science and technology by investigating the application of 
*science to technology and the impact of technology on science, e.g. Archimedes' screw to "Il- 
lustrate the principle of oits use in water irrigation or the use of springs and similar *propulsion 
devices in toys 

• 

• Experiment 4: Stretching Springs 
• 

• Hang a spring vertically using a clamp and stand. Add weights to the end of the spring, Mea- 
sure the length of the spring after each additional weight is added. Calculate the ""extension 
(Note: some initial weight will probably need to be added before the spring *stretches at all) 
Graph extension vs applied force. Interpret the graph, e.g. find the "Extension for unit applied 
force or applied force per unit extension. 



• •••••• investigate and understand the turning effect of a force. 

• Calculate the turning effect of a force = F x d where force (F) and distance (d) 
are ********* a j j-jgj^t angles. 

• Use the lever rule: Fldl = F2d2 • 

Discuss and brainstorm ideas about what causes things to turn about a pivot (fulcrum). 

• • • • • • Introduce the concept of the turning effect of a force (torque). 

• • Discuss the concept levers - to increase the turning effect, use a greater distance 
******* (e.g. a long crowbar) and/or a greater force (e.g. a heavier weight) 

• • Use the lever rule to calculate forces and distances at right angles to the applied 
********force (given three quantities and the fourth unknown). 
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• Experiment 5: The turning effect of a force 

• Kit available from the science Office 

•••••• Homework: Coursebook - The Turning Effect of a Force 

• Extension 



2.16 Changing mass 



So far we've assumed that the mass of the objects being considered is constant, which is not 
always true. Mass is conserved overall, but it can be useful to consider objects, such as rockets, 
which are losing or gaining mass. 

We can work out how to extend Newton's second law to this situation by considering a rocket 
two ways, as a single object of variable mass, and as two objects of fixed mass which are being 
pushed apart. 

We find that 



d{mv) 

F= — j — 
dt 

Force is the rate of change of a quantity, m v, which we call linear momentum. 



2.17 Newton's third law 

Newton's third law says that the force, F12 exerted by amass, m\, on a second mass, m 2 , is equal 
and opposite to the force F21 exerted by the second mass on the first. 

<j(m 2 v 2 ) _ p __ p djmivi) 
dt " dt 

if there are no external forces on the two bodies. We can add the two momenta together to get, 

d{m\Vi) d(m 2 v 2 ) d 
— — — + — — — = — (mivi + m 2 v 2 ) = 
dt dt dt 

so the total linear momentum is conserved. 

Ultimately this is consequence of space being homogeneous. 



2.18 Centre of mass 

Suppose two constant masses are subject to external forces, Fi, and F 2 
Then the total force on the system, F, is 
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F = Fi + F 2 = — (mivi + m 2 v 2 ) 
at 

because the internal forces cancel out. 

If the two masses are considered as one system, F should be the product of the total mass and , 
the average acceleration, which we expect to be related to some kind of average position. 

(mi + m 2 ) a = 4 (miVi + m 2 v 2 ) 
= ^2(miri + m 2 r 2 ) 

p. _ £_ m 1 r 1 + m 2 r 2 

u dt 2 rrii+mz 

The average acceleration is the second derivative of the average position, weighted by mass. 

This average position is called the centre of mass, and accelerates at the same rate as if it had the 
total mass of the system, and were subject to the total force. 

We can extend this to any number of masses under arbitrary external and internal forces. 

• we have n objects, mass mi, m 2 ... m n , total mass M 

• each mass m ; is at point rj, 

• each mass m\ is subject to an external force F ; 

• the internal force exerted by mass mj on mass m ; is F y 

then the position of the centre of mass, R, is 

R= Li mjXi 
M 

We can now take the second derivative of R 

mr = Zimn 
= LtPt + LjVij) 

But the sum of all the internal forces is zero, because Newton's third law makes them cancel in 
pairs. Thus, the second term in the above equation drops out and we are left with: 

MR=£F, =F 

i 

The centre of mass always moves like a body of the same total mass under the total external force, 
irrespective of the internal forces. 

<p> If a rigid body is initially at rest, it will remain at rest if and only if the sum of all the forces 
and the sum of all the torques acting on the body are zero. As an example, a mass balance with 
arms of differing length is shown in figure above. 

The balance beam is subject to three forces pointing upward or downward, the tension T in the 
string from which the beam is suspended and the weights Mig and M 2 g exerted on the beam 
by the two suspended masses. The parameter g is the local gravitational field and the balance 
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Centre of mass 



beam itself is assumed to have negligible mass. Taking upward as positive, the force condition 
for static equilibrium is 

T - Mi g - M 2 g = Zero net force 

Defining a counterclockwise torque to be positive, the torque balance computed about the pivot 
point in figure 10.7 is 

t = Migdi - M 2 gd 2 = Zero torque 
where d\ and d 2 are the lengths of the beam arms. 

The first of the above equations shows that the tension in the string must be 

T={M 1+ M 2 )g 

while the second shows that 

Mi _ d\ 

M 2 d 2 

Thus, the tension in the string is just equal to the weight of the masses attached to the balance 
beam, while the ratio of the two masses equals the inverse ratio of the associated beam arm 
lengths. 
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3 Rotational Dynamics 



There are two ways to multiply two vectors together, the dot product and the cross product. We 
have already studied the dot product of two vectors, which results in a scalar or single number. 

The cross product of two vectors results in a third vector, and is written symbolically as follows: 



Ax B 

The cross product of two vectors is defined to be perpendicular to the plane defined by these 
vectors. However, this doesn't tell us whether the resulting vector points upward out of the plane 
or downward. This ambiguity is resolved using the right-hand rule: 

1. Point the uncurled fingers of your right hand along the direction of the first vector A. 

2. Rotate your arm until you can curl your fingers in the direction of the second vector B. 

3. Your stretched out thumb now points in the direction of the cross product vector Ax B. 

The magnitude of the cross product is given by 

|AxB| = |A||B|sin0 

where | A| and | B| are the magnitudes of A and B, and 6 is the angle between these two vectors. 
Note that the magnitude of the cross product is zero when the vectors are parallel or anti-parallel, 
and maximum when they are perpendicular. This contrasts with the dot product, which is max- 
imum for parallel vectors and zero for perpendicular vectors. 

Notice that the cross product does not commute, i. e., the order of the vectors is important. In 
particular, it is easy to show using the right-hand rule that 

AxB=-BxA 

An alternate way to compute the cross product is most useful when the two vectors are expressed 
in terms of components, 



C=AxB= 



A x Ay A z 

B X By B z 



where the determinant is expanded as if all the components were numbers, giving 

C x = A y B z — A z B y 
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Rotational Dynamics 



are corresponding expressions for the other components. 

With the cross product we can also multiply three vectors together, in two different ways. 
We can take the dot product of a vector with a cross product, a triple scalar product, 

A-(BxC) = (AxB)-C 

The absolute value of this product is the volume of the parallelpiped defined by the three vectors, 
A, B, and C 

Alternately, we can take the cross product of a vector with a cross product, a triple vector product, 
which can be simplified to a combination of dot products. 

A x (B x C) = (A • C)B - (A • B)C 

This form is easier to do calculations with. 
The triple vector product is not associative. 

Ax(BxC)^(AxB)xC 
A nice as well as a useful way to denote the cross product is using the indicial notation 

C=AxB= e ijk AjB k ei 

where e l i k is the Levi-Civita alternating symbol and e; is either of the unit vectors i, j, k. (A good 
exercise to convince yourself would be to use this expression and see if you can get C = Ax B as 
defined before.) 

3.1 Torque 

Torque is the action of a force on a mass which induces it to revolve about some point, called the 
origin. It is defined as 

t = r x F 

where r is the position of the mass relative to the origin. 

Notice that the torque is zero in a number of circumstances. If the force points directly toward 
or away from the origin, the cross product is zero, resulting in zero torque, even though the force 
is non-zero. Likewise, if r=0, the torque is zero. Thus, a force acting at the origin produces no 
torque. Both of these limits make sense intuitively, since neither induces the mass to revolve 
around the origin. 
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3.2 Angular momentum 

The angular momentum of a mass relative to a point O is denned as 

L = r x p 

where p is the ordinary momentum of the mass. The angular momentum is zero if the motion of 
the object is directly towards or away from the origin, or if it is located at the origin. 

If we take the cross product of the position vector and Newton's second law, we obtain an equa- 
tion that relates torque and angular momentum: 

dp d , . dr 
rxF = rx — = — rxp — -xp 

dt dt y ' dt 

Since the cross product of parallel vectors is zero, this simplifies to 

dL 

T= — 

dt 

This is the rotational version of Newton's second law. 

For both torque and angular momentum the location of the origin is arbitrary, and is generally 
chosen for maximum convenience. However, it is necessary to choose the same origin for both 
the torque and the angular momentum. 

For the case of a central force, i. e., one which acts along the line of centers between two objects 
(such as gravity), there often exists a particularly convenient choice of origin. If the origin is 
placed at the center of the sun (which is assumed not to move under the influence of the planet's 
gravity), then the torque exerted on the planet by the sun's gravity is zero, which means that the 
angular momentum of the planet about the center of the sun is constant in time. No other choice 
of origin would yield this convenient result. 

We already know about two fundamental conservation laws -- those of energy and linear mo- 
mentum. We believe that angular momentum is similarly conserved in isolated systems. In other 
words, particles can exchange angular momentum between themselves, but the vector sum of 
the angular momentum of all the particles in a system isolated from outside influences must 
remain constant. 

In the modern view, conservation of angular momentum is a consequence of the isotropy of 
space -- i. e., the properties of space don't depend on direction. This is in direct analogy with 
conservation of ordinary momentum, which we recall is a consequence of the homogeneity of 
space. 

3.3 Angular velocity and centrifugal force 

If an object is rotating about an axis, n, n being a unit vector, at frequency to we say it has angular 
velocity to. Despite the name, this is not the rate of change of a angle, nor even of a vector. 
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If a constant vector r is rotating with angular velocity to about a fixed point then 

f = o) x r 

This says the acceleration is always at right angles to both the velocity and the axis of rotation. 

When the axis is changing a) can be defined as the vector which makes this true. 

Note that on the left hand side of this equation r is a vector in a fixed coordinate system with 
variable components but on the right hand side its components are given in a moving coordinate 
system, where they are fixed. 

We can distinguish them more clearly by using subscripts, 'r' for rotating and T for fixed, then 
extend this to arbitrary vectors 

g/ = gr+wxg / 

for any vector, g. 

Using this, we can write Newton's second law in the rotating frame. 

V f = mvf 

= mv r + mwxvf 

= mr r + mft)x(2r r + d)xr|) 

or, rearranging 

ma r = F-2m(w x v r ) - mo) x(wxr) 

The mass behaves as if there were two additional forces acting on it. The first term, a> x v is called 
the Coriolis force. The second term is recognizable as the familiar centrifugal force. 

The kinetic energy and angular momentum of the dumbbell may be split into two parts, one 
having to do with the motion of the center of mass of the dumbbell, the other having to do with 
the motion of the dumbbell relative to its center of mass. 

To do this we first split the position vectors into two parts. The centre of mass is at. 

= MiT] + M 2 r 2 
Mi + M 2 

so we can define new position vectors, giving the position of the masses relative to the centre of 
mass, as shown. 

r*=ri-Ri r*=ri-Ri 

The total kinetic energy is 
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K = \MiVf + \MiVl 

= iMilV-v*! 2 + ±M 2 |V-v*| 2 

= ±(Mi + M 2 )V 2 + \Mi_vf- + \M 2 vf 

- K ext + Ki nt 

which is the sum of the kinetic energy the dumbbell would have if both masses were concen- 
trated at the center of mass, the translational kinetic energy and the kinetic energy it would have 
if it were observed from a reference frame in which the center of mass is stationary, the rotational 
kinetic energy. 

The total angular momentum can be similarly split up 

= (Mi + M 2 )R x V + (Mir * x v* + M 2 r* x v* ) 

into the sum of the angular momentum the system would have if all the mass were concentrated 
at the center of mass, the orbital angular momentum, and the angular momentum of motion 
about the center of mass, the spin angular momentum. 

We can therefore assume the centre of mass to be fixed. 

Since a> is enough to describe the dumbbells motion, it should be enough to determine the 
angular momentum and internal kinetic energy. We will try writing both of these in terms of a) 

First we use two results from earlier 



L = r x v and v = wxr 



to write the angular momentum in terms of the angular velocity 

L = Mir* xv* + M 2 r^ x v^ 

Mir*x(a>xr*) + M 2 r*x(wxr*) 

= MiMr*.r*)-r*(r*.a/)) + M 2 (a>(r* -r*)-r* (r* ■«>)) 

Mi [todf - r* (r* • w)) + M 2 {bid\ - r* (r* • o>)) 

(M\d\ + M 2 d\)oi - (Mir*(r*-w) + M 2 r*(r*-w)) 

The first term in the angular momentum is proportional to the angular velocity, as might be 
expected, but the second term is not. 

What this means becomes clearer if we look at the components of L For notational convenience 
we'll write 



A = (xi.yi.zi) r* 2 = {x2,y 2 ,z 2 ) 
These six numbers are constants, reflecting the geometry of the dumbbell. 
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L x = {M 1 d 2 + M 2 d 2 2 )a) x -(M l x 2 + M 2 x 2 )a) x 

- (Mi xi yi + M 2 x 2 y 2 ) 0) y ~ {Mi x L Z\ + M 2 x 2 z 2 ) a) z 

L y = {M\d\ + M 2 dl)o) y -{Mix\ + M 2 xl)(i)y 

-{MiXiyi + M 2 x 2 y 2 )(D x -{M l yiZi + M 2 y 2 z 2 )(D z 

L z = {M 1 d 2 + M 2 d 2 )a) z -(MiX 2 + M 2 x 2 )a) z 

-(M1X1Z1 + M 2 x 2 z 2 )(D x - {M 1 y 1 z l + M 2 y 2 z 2 )(o y 

This, we recognise as being a matrix multiplication. 



(LA 




'I XX 


Ixy 


h 




lyx 


hy 


\L ZJ 




Jzx 


Izy 



y z 









Q)y 


J 





where 






^xy 


lyx 


hy 


Jzx 


Izy 



Ixz\ 

lyz 

IzzJ 



{M x {d\-x\) + M 2 {d 2 -x. 
-(M 1 x 1 yi+M 2 x 2 y 2 ) 
-(M1X1Z1 +M 2 x 2 z 2 ) 



-{M\Xiyi + M 2 x 2 y 2 ) 
M x {d\-y\) + M 2 {d\ 
-{M x y x z y + M 2 y 2 z 2 ) 



y 2 2 ) 



The nine coefficients of the matrix I are called moments of inertia. 
By choosing our axis carefully we can make this matrix diagonal. E.g. if 



-(M!XiZi + M 2 x 2 z 2 ) 

-(Miyizi + M 2 y 2 z 2 ) 
M x {d\-z\) + M 2 {d\-z 



r* = (di,0,0) r* = (-d 2 ,0,0) 

then 



'Ixx 


Ixy 


Ixz\ 


fo 





\ 


lyx 


hy 


Iyz\ = 





M Y d\ + M 2 d\ 





\Izx 


Izy 


Izz) 


lo 





M\d\ + M 2 d\ t 



Because the dumbbell is aligned along the x-axis, rotating it around that axis has no effect. 
The relationship between the kinetic energy, T, and oj quickly follows. 



2T = M\ v* 2 + 

= M!V*-(o>xr*) + MiV* 2 -{(D*x* 2 ) 
= M^-fr* xv*) + M 2 (i)-(r*xv*) 



M 2 v* 2 



On the right hand side we immediately recognise the definition of angular momentum. 



T= -a) L 

2 



Substituting for L gives 



T = -tola) 

2 
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Further reading 



Using the definition 



a) = wn n • n = 1 



this reduces to 



1 7 
2 

where the moment of inertia around the axis 'n is 



/ = nln 



a constant. 

If the dumbbell is aligned along the x-axis as before we get 

T=^{M l df + M 2 dZ)(oj 2 y + ct) 2 z ) 

These equations of rotational dynamics are similar to those for linear dynamics, except that I is 
a matrix rather than a scalar. 



3.4 Further reading 

• Modern Physics/The Uneven Dumbbell 



3.5 Summation convention 



If we label the axes as 1,2, and 3 we can write the dot product as a sum 

3 

U V = £ UiVi 
i=l 

If we number the elements of a matrix similarly, 



A = 



'A u An A l3 ' 
A21 A22 A23 
,A 3i A 32 A 33 ) 



B 



'Bn Bn #13' 
B21 B22 B23 

Al B32 #33; 



we can write similar expressions for matrix multiplications 



1 http://en.wikibooks.org/wiki/Modern%20Physics%2FThe%20Uneven%20Dumbbell 
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3 



3 



{Au)i = Y J A ij u j (AB)/jfc 



Notice that in each case we are summing over the repeated index. Since this is so common, it is 
now conventional to omit the summation sign. 

Instead we simply write 

u v = tyVi (Au)j = AijUj (AB)j fc = AijBjk 
We can then also number the unit vectors, e and write 

U = U{%i 

which can be convenient in a rotating coordinate system. 

3.6 Kronecker delta 

The Kronecker delta is 




This is the standard way of writing the identity matrix. 



3.7 Levi-Civita (Alternating) symbol 



Another useful quantity can be defined by 



f 1 [i.j.k) 
eijk = \-l d,j,k) 




(1,2,3) or (2,3,1) or (3,1,2) 
(2,1,3) or (3,2,1) or (1,3,2) 
otherwise 



With this definition it turns out that 



uxv = e ijk eiU]V k 



and 



e ijk e ipq — djpfikq <5 jqSkp 



26 



Levi-Civita (Alternating) symbol 



This will let us write many formulae more compactly. 

The uneven dumbbell consisted of just two particles. These results can be extended to cover 
systems of many particles, and continuous media. 

Suppose we have AT particles, masses mi to m n, with total mass M. Then the centre of mass is 

_ ^Ln m n r n 
~ M 

Note that the summation convention only applies to numbers indexing axes, not to n which 
indexes particles. 

We again define 

r^ = r„-R 

The kinetic energy splits into 

T = l -MV 2 + Wm n v* n 2 
z z n 

and the angular momentum into 

L = MR xV+^ m n r* n x v* 

n 

It is not useful to go onto moments of inertia unless the system is approximately rigid but this is 
still a useful split, letting us separate the overall motion of the system from the internal motions 
of its part. 

If the set of particles in the previous chapter form a rigid body, rotating with angular velocity oj 
about its centre of mass, then the results concerning the moment of inertia from the penultimate 
chapter can be extended. 

We get 

hj = Y, m n( r l s ij - r mr n j) 

n 

where ( r n \, r n 2, r ^3) is the position of the n th mass. 
In the limit of a continuous body this becomes 

Iij = f v p(.r){r 2 dij-rirj)dV 

where p is the density. 

Either way we get, splitting L into orbital and internal angular momentum, 
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Li = Me ijk R j V k + I ij o) j 
and, splitting T into rotational and translational kinetic energy, 

T=^MV i V i + ^(o i I ij (Oj 
It is always possible to make I a diagonal matrix, by a suitable choice of axis. 



4 Mass Moments Of Inertia Of Common 
Geometric Shapes 



The moments of inertia of simple shapes of uniform density are well known. 

4.0.1 Spherical shell 

mass M, radius a 



4.0.2 Solid ball 

mass M, radius a 



2 , 

I xx = I yy = hz = -Ma z 



I xx = Iyy = hz = 7 Ma 



4.0.3 Thin rod 

mass M, length a, orientated along z-axis 

1 



Ixx = Iyy=— Ma 2 I zz = 



4.0.4 Disc 

mass M, radius a, in x-y plane 



1 , 1 2 

I xx = I yy = -Ma z I zz = -Ma z 



4.0.5 Cylinder 

mass M, radius a, length h orientated along z-axis 



Mass Moments Of Inertia Of Common Geometric Shapes 




) 



1 



Ma 2 



4.0.6 Thin rectangular plate 



mass M, side length a parallel to x-axis, side length b parallel to y-axis 




) 



4.1 further reading 

• Statics /Moment of Inertia (contents) 1 

• Statics /Geometric Properties of Solids 2 



1 http://en.wikibooks. org/w i k i /Stat i c s %2FMoment%2 0of%2 Inert i a%2 0%28 content s%2 9 

2 http: //en. wikibooks.org/wiki /Stat i cs%2FGeometric%20Propert ie s%20of%20Sol ids 
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5 Newton's Laws: A second look 



So far we've tacitly assumed we can just calculate force as a function of position, set up the ODE 

raf = F(r) 

and start solving. 

It is not always so simple. 

Often, we have to deal with motion under constraint; a bead sliding on a wire, a ball rolling 
without slipping, a weight dangling from a string. 

There has to be some force keeping the bead on the wire, but we don't know what it is in advance, 
only what it does. This isn't enough information for us to write down the ODE. 

We need a way of solving the problem without knowing the forces in advance. 

How easy this is depends on the type of constraint. 

• If the constraint is an inequality, as with the weight on the string, there is no straightforward 
analytical method. 

• If the constraint can be written as a set of differential equations, and those equations can't be 
integrated in advance, there is an analytical method, but it is beyond the scope of this book. A 
ball rolling without slipping falls into this category. 

• If the constraint can be written as a set of algebraic equations, and frictional forces are negli- 
gible, there is a straightforward method that solves the problem. 

5.1 Generalised coordinates 

Suppose we have a system of n particles satisfying k constraints of the form 

/fc(ri,---,r„,f) = 

then we can use the constraints to eliminate k of the 3 n coordinates of the particles, giving us a 
new set of 3 n- k independent generalised coordinates; q\, qt,... q3 n-'k </sub>. 

Unlike the components of the position vectors, these new coordinates will not all be lengths, and 
will not typically form vectors. They may often be angles. 

We now need to work out what Newton's laws will look like in the generalised coordinates. 
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5.1.1 Derivation 

The first step is to eliminate the forces of constraint. 

We will need to consider a virtual displacement. This is an infinitesimal displacement made, 
while holding the forces and constraints constant. It is not the same as the infinitesimal displace- 
ment made during an infinitesimal time, since the forces and constraints may change during 
that time. 

We write the total force on particle i as 

F, =Ff +F- 

the sum of the externally applied forces and the forces of constraint. 
Newton's second law states 

We take the dot product of this with the virtual displacement of particle i and sum over all 
particles. 

£(F? + F<-p,Hr, = 

i 

We now assume that the forces of constraint are perpendicular to the virtual displacement. This 
assumption is generally true in the absence of friction; e.g, the force of constraint that keeps a 
ball on a surface is normal to the surface. 

This assumption is called D'Alembert's principle. Using it we can eliminate the forces of con- 
straint from the problem, giving 

L(F?-*/)-ffr / = 

i 

or 

; i 

The left hand side of this equation is called the virtual work. 
Now we must change to the generalised co-ordinate system. 
We write 

*i = r;(<7i,<72>" - ,qzn-k> i) 

Using the chain rule gives 
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Generalised coordinates 



and 



j dqj 1 dt 



Note that this implies 



dqj dqj 

The virtual work is, dropping the superscript, 



(2) 



LjQjSqj 

where the Qj are the components of the generalised force. 
We now manipulate the right hand side of (1) into a form comparable with this last equation 
The right hand side term is 



LiPi-Sri = Lirriiii-Sri 
The terms in the coefficient of qj can be rearranged 



= LtjmiTf^jSq] 



on substituting in equation (2) from above 
On taking a close look at this last equation, we see a resemblance to the total kinetic energy, 

i 

We now further rearrange to get an expression explicitly involving T. 



- <L 

dt 



d ( d j] _ d t 
dt \dcjj j dqj 



Putting this last expression into (1) along with the generalised force give 



33 



Newton's Laws: A second look 



E 



d 


[—) 


dT 1 


~dt 


{dqj) 





8qj = 



Since the 6 qj, unlike the 8 r are independent, this last equation can only be true if all the 
coefficients vanish. 

That is we must have 

d dT dT 

= Qi (3) 

dtdqj dqj 

These are the equations of motion for the system, in a general set of coordinates for which all 
constraints are automatically satisfied. 

For example, suppose we have a cylinder, mass m, radius a, rolling without slipping on a flat 
plane. 

The kinetic energy of the cylinder is 

T = - mx 2 + -ma 2 d 2 
2 8 

using the results from Rigid Bodies 1 , where x is the axis in the plane perpendicular to the axis 
of the cylinder, and 6 is the angle of rotation. 

Rolling without slipping implies 

x= ad 

so we get 



5 2 5 
T=-mx -mx = Q x 
8 4 v 

The cylinder has the same kinetic energy as if its mass were 20% greater. If there is no torque 
on the cylinder then Q x = F x , and the cylinder behaves in every respect as though it were a 20% 
larger point mass. 

To use (3) more generally, we need an expression for the Qj 
Suppose, as is often the case, that 

d 

then, by definition 



1 Chapter 3 . 7 on page 2 7 



34 



Generalised coordinates 



LjQjSqj = 




so, equating coefficients, the generalised force is 



Qj = 



dV 



Putting this generalised force into (3) gives 



d d{T-V) 



d(T-V) 
dqj 



= 



dt ddj 



since V has been assumed independent of the velocities. 

In fact, this last equation will still be true for some velocity dependent forces, most notably mag- 
netism, for a suitable definition of V, but we won't prove this here. 

We call the T-Vthe Lagrangian, I, and write 



We call these equations Lagrange's equations. They are useful whenever Cartesian co-ordinates 
are inconvenient, including motion under constraint. 

5.1.2 Example 

Suppose we have two identical point masses, m, connected by a string, length a. The string is 
threaded through an hole in a flat table so that the upper mass is moving in a horizontal plane 
without friction, and the lower mass is always vertically below the hole. The distance of the upper 
mass from the hole is r. 

The position of the mass on the table is best described using polar coordinates, ( r,8). Its kinetic 
energy is then 



d dL dL 




dtddj dqj 



-m[r 2 + r 2 2 ) 



The velocity of the lower mass is d(a-r)/dt= -drldt, so the total kinetic energy is 



T= -m{2f 2 + r 2 9 2 ) 



The potential energy is 
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V = mg[a-r) 



where g is the gravitational acceleration. 



This means 



L = -m[2r 2 + r 2 6 2 ) + mg{r - a) 



and the equations of motion are 



d dL 



dL 
dO 
dL 
dr 



mf t {r 2 6) 



= 




= 2mr-mr6 2 



+ mg = 



dt dr 



The first of these equations says that the angular momentum is constant, as expected since there 
is no torque on the particles. If we call this constant angular momentum I then we can write 



mr 2 

and the second equation of motion becomes 



then the lower ball will be pulled out of the hole, at which point these equations of motions cease 
to apply. They only hold when 0< r< a, a constraint which is not easily tamable. 

Notice, we have not needed to calculate the tension in the string, which is the force of constraint 
in this problem. 

5.2 Reformulating Newton 

In the last chapter we saw how to reformulate Newton's laws as a set of second order ordinary 
differential equations using arbitrary coordinates: 



r = 




2m 2 r 3 2 



Clearly, if initially 



I 2 > gm 2 r 3 , 



d dL dL 



dt dcji dqi 
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Derivation 



It is often far easier to solve a problem by starting with the Lagrangian than by explicitly using 
Newton's laws. The Lagrangian is also more useful for theoretical analysis. 

However, the Lagrangian approach is not the only useful way to reformulate Newton. There is a 
second, related, approach which is also highly useful for analysis; namely, writing the equations 
of motion as a set of first order ordinary differential equations in a particularly natural way. 

To do this we will use a standard technique, familiar from your calculus studies. 



5.3 Derivation 

Consider the function 



H{q,p, t) = qiPi-Uq, q, t) 

where the q ; are the generalised coordinates, the p ; are a set of new variables whose meaning 
we will soon see, H is to be a function of the p's and q's alone, and the summation convention 
is being used. 

Then we have 

dL J dL J dL J 
dH = qidpi + pidqi - ^—dqi - - — dqi - —at (1) 
dqt dqi dt 

but, since H is a function of the p's and q's alone we can also write 

dH , dH , dH , 

dH= dqt + dpi + — dt (2) 

dqi dpi dt 

For these two equations to both be true the coefficients of the differentials must be equal. 

Equation (2) does not contain a term in dqi, so the coefficient of that term in equation (1) must 
be zero. I.e, 

dL 

dqi 

which gives us a definition of the p ,-. Using this in Lagrange's equations gives 

dL 

dqi 

and equation (1) simplifies to 

dL 

dH= qidpi- pidqi - — dt 
dt 
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A comparison of coefficients with equation (2) now gives the desired set of first order equations 
for the motion, 

dH dH 
dpi dm 

These are called Hamilton's equations and H is called the Hamiltonian. 



5.4 Physical meaning of H and p 

To see what H and the p ;'s actually are, let's consider a few typical cases. 

First, let's look at a free particle, one subject to no forces. This will let us see what the p mean, 
physically. 

If we use Cartesian coordinates for the free particle, we have 

L=T=^m{x 2 + f + z 2 ) 
The p , 's are the differential of this with respect to the velocities. 

dL dL dL 

v*=j-± = mk p y = ry = mf Pz = di = mi 

These are the components of the momentum vector. 
If we use cylindrical coordinates, we have 

L=T=^m{f 2 + r 2 9 2 + z 2 ) 

and 

dL dL 2 . dL 

p r = — = mr pa - — = mr 9 p z = — = mz 
df pw 39 dz 

This time, p z is the component of momentum is the z direction, p r is the radial momentum, 
and po is the angular momentum, which we've previously seen is the equivalent of momentum 
for rotation. 

Since, in these familiar cases, the p ;'s are momenta, we generalise and call the p/s conjugate 
momenta. 

Note that from Hamilton's equations, we see that if H is independent of some coordinate, q, 
then 

dH 
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Deeper meaning 



so the momentum conjugate to q is conserved. This connection between conservation laws and 
coordinate independence lies at the heart of much physics. 

Secondly, to see the physical meaning of H, we'll consider a particle moving in a potential V, 
described using Cartesian coordinates. 

This time the Lagrangian is 

L=T-V = ^m[x 2 + y 2 + z 2 )-V{x,y,z) 

We find that the p ,'s are 

dL dL dL 

Px = ^r = mx p y = — = my p z = — = mz 
ox dy oz 

so, again, they are the momenta 
Now lets calculate H. 



H = xp x + ypy + yp y - L 

m[x 2 + y 2 + z 2 ) - \m[x 2 + y 2 + z 2 ) + V 

= \ m [x 2 + y 2 + z 2 ) + V 

T + V 

so H is the total energy, written as a function of the position and momentum. 

This is not always true, but it will be true for all common systems which conserve energy. 

In general, if we know how the energy of a systems depends on speeds and positions, we know 
everything about the system. 

Physicists will often start work on a problem by writing down an expression for H, or L, and 
never bother calculating the actual forces. 



5.5 Deeper meaning 

If we compare Hamilton's equations, 

dH dH 

dpi dqi 

with the equations of geometrical optics. 

. _ do ^ dct) 
1 dki 1 dxi 

we see that both sets of equations have the same form. 
If we made a substitution of the form, 
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H=a<x> p = ak 



for any a, the two sets of equations would become identical. 

Now, geometrical optics is an approximation to the full solution of any wave equation, valid 
when the wavelength is small. 

Since classical mechanics is just like geometrical optics, classical mechanics could also be small 
wavelength approximation to a more accurate wave theory. 

This is not a proof, could be is not the same as is, but it does show that classical mechanics is 
compatible with wave theories, contrary to intuition. 

We will see later that this is not a coincidence, classical mechanics really is a small wavelength 
approximation to a more accurate wave theory, namely quantum mechanics, and that, in that 
theory, the energy of matter waves is proportional to frequency and the momentum to wave 
number, just as the substitution above requires. 

Hamilton's equations can also be used to construct a mathematical gadget, called a Poisson 
bracket, which lets us write the equations of classical mechanics in a way which makes the 
connections with quantum mechanics even clearer, but an explanation of how would be beyond 
the scope of this book. 



<! — A tease, but it should encourage the reader to look deeper — > 



Since classical mechanics is just like geometrical optics, there should be an analog to Fermat's 
principle, that light takes the quickest path, and there is. 

Where light takes the quickest path from a to b, matter takes the path of least action, where the 
action is defined as 



It can be proved that this integral takes extremal values if and only if the system obeys Lagrange's 
equations. Saying matter takes the path of least action is equivalent to saying it obeys Lagrange's 
equations, or Newton's laws. 

Furthermore, it turns out that in general relativity, the action is proportional to the time taken, 
so matter is also really taking the quickest path between two points. 

We've seen that an Hamiltonian of the form 

H= ^ + V(r) 
2m 

describes motion under a conservative force, but this is not the most general possible Hamilto- 
nian. 

Many waves are described by another simple form, in the geometrical optics limit, 




40 



Deeper meaning 



H = /(|p|) 

What happens if we consider other forms for the Hamiltonian? 
Suppose we have 

PP 

H= — +p-A(r) 

The dot product term is the simplest scalar we can add to the kinetic energy that's dependent on 
momentum. We'll see that it acts like a potential energy. 

Using Hamilton's equations we can immediately write down equations of motion. 

r- - M - Ei+A- 

- _dM. — _ £A/ 
Pi ~ dxi ~ PJ dxt 

Note that we are using the summation convention 2 here. 

Now that we've got the equations of motion, we need to work out what they mean. 

The first thing we notice is that the momentum is no longer ra'v . There is an additional con- 
tribution from the potential field 'A. We can think of this as a form of potential momentum, 
analogous to potential energy. 

The force on the particle is 

d 

Fi = mxi = pi + m—Aiix) 
at 

Using the chain rule, and substituting for dpi dt with the second equation of motion, we get 

OAj . dAi 
F i = - pj — + mxj — 

Using the first equation, this becomes 



F{ = mxj 



dAi dA i\ dA i 



, dxj dxi j 1 dx, 



The term in brackets is recognizable as being the type of thing we see in cross products. With a 
little manipulation, using the Kronecker delta and alternating symbol, we can write 



2 http://en.wikibooks.org/wiki/. .%2FIndex%20Notation 
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Ft = tjiStiSjk-ettSjfitfe + mAj^ 

dAi . dAj 

Xj£ n ije n iicg^ + mA j~d^- 

-eij„Xj(V xA) n + mAj-g^- 

-EijnXj (V x A) n + \-fcAjAj 



so, writing the last line as a vector equation, 



1 



2 



F = -mvx (V x A) + -mVA 
2 

The force from this potential has a component perpendicular to its curl and to the velocity, and 
another component which is the gradient of a scalar. 

If we add a carefully chosen potential energy to the Hamiltonian we can cancel out this second 
term. 

P 2 1 7 

H = — +p-A{r) + -mA z F = -mvx(VxA) 

2m 2 
This Hamiltonian can be simplified to 

H= — fp + mA) 2 
2m ^ ' 

where the term in bracket is m'v , written as a function of momentum. 

This simple modification to the Hamiltonian gives us a force perpendicular to velocity, like mag- 
netism. Because the force and velocity are perpendicular, the work done by the force is always 
zero. 

We can use potential fields to describe velocity dependent forces, provided that they do no work. 

The coefficient of A in these expression is arbitrary. Changing it amount merely to measuring A 
in different units so we can equally well write 

H= (p + aA) F=-avx(VxA) 

Having a= m means A is measured in units of velocity, which may sometimes be convenient, 
but we can use any other constant value of a that suits our purpose. When we come to study 
relativity, we'll use a=-l. 

Notice that the force depends only on the curl of A, not on A itself. This means we can add any 
function with zero curl to A without changing anything, just as we can add a constant to the 
potential energy. 

Furthermore, it is a standard result of vector calculus that any vector field is the sum of two 
components, one with zero curl, the other with zero divergence. Since the zero curl component 
does not affect the total force, we can require it to be zero; i.e. we can require A to have zero 
divergence. 
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V-A = 

This is called a gauge condition. For the moment, it can be considered to define especially 
natural values for the integration constants. 

Mach's principle 3 

Machs principle is the name given by Einstein to a vague hypothesis first supported by the physi- 
cist and philosopher Ernst Mach. 



3 http://en.wikipedia.org/wiki/Mach%27s%20principle 
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6 Harmonic Oscillators 



6.1 Energy Analysis 



Figure 11.2: Potential, kinetic, and total energy of a one-dimensional harmonic oscillator plot- 
ted as a function of spring displacement. 

For a spring, Hooke's law says the total force is proportional to the displacement, and in the 
opposite direction. 



Since this is independent of velocity, it is a conservative force. We can integrate to find the po- 
tential energy of the mass-spring system, 

V{x) = \kx 2 

Since a potential energy exists, the total energy 



is conserved, i. e., is constant in time. 

We can now use energy conservation to determine the velocity in terms of the position 



We could integrate this to determine the position as a function of time, but we can deduce quite 
a bit from this equation as it is. 

It is fairly evident how the mass moves. From Hooke's law, the mass is always accelerating toward 
the equilibrium position, so we know which sign of the square root to take. 

The velocity is zero when 



If x were larger than this the velocity would have to be imaginary clearly impossible, so the mass 
must be confined between these values. We can call them the turning points. 



F = -kx 



1,2 1 ? 
E= - kx + - mx 
2 2 





45 



Harmonic Oscillators 



If the mass is moving to the left, it slows down as it approaches the left turning point. It stops 
when it reaches this point and begins to move to the right. It accelerates until it passes the equi- 
librium position and then begins to decelerate, stopping at the right turning point, accelerating 
toward the left, etc. The mass thus oscillates between the left and right turning points. 

How does the period of the oscillation depend on the total energy of the system? We can get a 
general idea without needing to solve the differential equation. 

There are only two parameters the period, T, could depend on; the mass, m, and the spring 
constant, k. 

We know T is measured in seconds, and m in kilograms. For the units in Hooke's law to match, 
/cmust be measured in N-m" 1 , or equivalently, in kg-s" 2 . 

We immediately see that the only way to combine m and k to get something measured in sec- 
onds is to divide, cancelling out the kg's. 

Therefore Toe V mlk 

We have established the general way the period depends on the parameters of the problem, with- 
out needing to use calculus. 

This technique is called dimensional analysis, and has wide application. E.g., if we couldn't 
calculate the proportionality constant exactly, using calculus, we'd be able to deduce by doing 
one experiment. Without dimensional analysis, if calculus failed us we'd have to do scores of 
experiments, each for different combinations of m and k. 

Fortunately, the proportionality constants are typically small numbers, like 3\/2 or 0T/2. 

6.2 Analysis Using Newton's Laws 

The acceleration of the mass at any time is given by Newton's second law 

_ d 2 x _ F _ kx 
dt 2 m m 

An equation of this type is known as a differential equation since it involves a derivative of the 
dependent variable . Equations of this type are generally more difficult to solve than algebraic 
equations, as there are no universal techniques for solving all forms of such equations. In fact, it 
is fair to say that the solutions of most differential equations were originally obtained by guess- 
ing] 

There are systematic ways of solving simple differential equations, such as this one, but for now 
we will use our knowledge of the physical problem to make an intelligent guess. 

We know that the mass oscillates back and forth with a period that is independent of the am- 
plitude of the oscillation. A function which might fill the bill is the sine function. Let us try 
substituting, 

x = Asincjt 
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where a) is a constant, into this equation. 
We get 

2 k 
-to Asmcjt = -A — sin&if 

m 

Notice that the sine function cancels out, leaving us with a) 2 = kl m. The guess thus works if we 
set 

(o = \ — 
V m 

This constant is the angular oscillation frequency for the oscillator, from which we infer the pe- 
riod of oscillation to be 

fm 

T - 2n U 

This agrees with the result of the dimensional analysis. Because this doesn't depend on A, we 
can see that the period is independent of amplitude. 

It is easy to show that the cosine function is equally valid as a solution, 

x = B coswf 

for the same a>. 

In fact, the most general possible solution is just a combination of these two, i. e. 

x = Asincjt + B cosa)t 

The values of A and B depend on the position and velocity of the mass at time f=0. 

If we wiggle the left end of the spring in the above diagram by amount d= dgsm a) F</sub>t, 
rather than rigidly fixing it, we have a forced harmonic oscillator. 

The constant do is the amplitude of the imposed wiggling motion. The forcing frequency a> p 
is not necessarily equal to the natural or resonant frequency of the mass-spring system. Very 
different behavior occurs depending on whether it is less than, equal to, or greater than a). 

Given the above wiggling, the force of the spring on the mass becomes 

F = - k(x-d) = -k{x- dosincjpt) 

since the length of the spring is the difference between the positions of the left and right ends. 
Proceeding as for the unforced mass-spring system, we arrive at the differential equation 
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d x k k 
~ 77 + — x = — dosmtopt 
dt 2 M M 

The solution to this equation turns out to be the sum of a forced part in which xoc sin oj F<lsub>t, 
and a free part which is the same as the solution to the unforced equation. We are primarily 
interested in the forced part of the solution, so let us set x= do sin to F</sub>tand substitute this 
into the equation of motion, giving: 



k k 
-oy^xo sin cop t + — xosin<yp t= — dosino)pt 

The sine factor cancels leaving us with an algebraic equation for , the amplitude of the oscillatory 
motion of the mass. 

Solving for the ratio of the oscillation amplitude of the mass to the amplitude of the wiggling 
motion, we find 



xp_ _ 1 _ 1 

d ~ l-Mco 2 F /k ~ l-(o 2 F /o) 2 

where we have recognized that kl M=o> 2 , the square of the frequency of the free oscillation. 

Notice that if a) p<a>, then the motion of the mass is in phase with the wiggling motion and the 
amplitude of the mass oscillation is greater than the amplitude of the wiggling. As the forcing 
frequency approaches the natural frequency of the oscillator, the response of the mass grows in 
amplitude. 

When the forcing is at the resonant frequency, the response is technically infinite, though practi- 
cal limits on the amplitude of the oscillation will intervene in this case -- for instance, the spring 
cannot stretch or shrink an infinite amount. In many cases friction will act to limit the response 
of the mass to forcing near the resonant frequency. 

When the forcing frequency is greater than the natural frequency, the mass actually moves in the 
opposite direction of the wiggling motion -- i. e., the response is out of phase with the forcing. 
The amplitude of the response decreases as the forcing frequency increases above the resonant 
frequency. 

Forced and free harmonic oscillators form an important part of many physical systems. For in- 
stance, any elastic material body such as a bridge or an airplane wing has harmonic oscillatory 
modes. A common engineering problem is to ensure that such modes are damped by friction or 
some other physical mechanism when there is a possibility of excitation of these modes by nat- 
urally occurring processes. A number of disasters can be traced to a failure to properly account 
for oscillatory forcing in engineered structures. 

We often encounter systems which contain multiple harmonic oscillators, such as this: 
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Figure 1 

two identical masses, m, the first attached to a wall by a spring with constant k, and the second 
attached to the first by another, identical spring. 

If the springs weren't linked they'd both vibrate at the same frequency, co=^/{ kl m). Linking the 
springs changes this. 

To find out how the linked system behaves, we will start with the Lagrangian, using the displace- 
ments of the masses, X\ and x 2 , as our coordinates. 

A moment's inspection of the system shows 

T=^m{xj + xl) V=^k[x 2 + {x 2 -x 1 ) 2 ) 

so, using a) 2 = kl m, 

L = - m [x\ + x 2 ) - - ma> 2 [x\ + {x 2 - X\) 2 ) 
The equations of motion immediately follow. 
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xi = -a) 2 {2x l -x 2 ) ( ^ 

X 2 - -CD 2 {X 2 -Xl) 

To solve these equations we try a solution in trig functions 

Xi=AisinQ£ x 2 = A 2 sin£lt 

Substituting this into (1) gives 

-CL 2 Ai + co 2 {2A l -A 2 ) = 
-Q 2 A 2 + w 2 {A 2 -A l ) = 

We would get the same equations from any trig function solution of the same frequency 
Gathering the coefficients of A\ and A 2 together lets rewrite the last equation as 




We can only solve this equation if the determinant of the matrix is zero. 

n 4 -3o) 2 D 2 + a) i = 

The solutions are 

n 2 + = 3 -^a> 2 

± 2 

so the combined system has two natural frequencies, one lower and one higher than the natural 
frequency of the individual springs. This is typical. 

We can also calculate the ratio of A\ and A 2 from (2). Dividing by A 2 gives 

M =1 _«± 

A 2 0) 2 

• For the lower frequency, Q_, which is less than oj, the ratio is positive, so the two masses move 
in the same direction with different amplitudes. They are said to be in phase. 

• For the higher frequency, D.+, which is more than a>, the ratio is negative, so the two masses 
move in the opposite direction with different amplitudes. They are said to be in antiphase. 

This behaviour is typical when pairs of harmonic oscillators are coupled. 

The same approach can also be used for systems with more than two particles. 
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6.3 Introduction 

In principle, we could use the methods described so far to predict the behavior of matter by 
simply keeping track of each atom. 

In practice, this is not a useful approach. Instead, we treat matter as a continuum. 

In this section we will see how this is done, and how it leads us back to waves. 

For an example, we will consider a set of N+l identical springs and masses, arranged as in the 
previous section, with spring attached to the wall, and mass N free. 

We are interested in what happens for large N, the continuum limit. 

The system has three other parameters; the spring constant k, the particle mass Am, and the 
spring rest spacing Aa, where variable names have been chosen for future convenience. These 
can be combined with N to give a similar set of parameters for the system. 

Suppose all the masses are displaced by d from rest, with total displacement D= dN of the 
system end, then the total potential energy is kNd 2 /2= kD 2 / (2 N), so 

• The system's spring constant is K~ k/N 

• The system mass is m= AmN 

• The system rest length is o= \aN 

If, as we increase N we change the other three parameters to keep K, m, and a constant then in 
the large N limit this discrete system will look like a spring with mass continuously distributed 
along its length. 

For coordinates, we will use the displacements of each mass, x „. For large N the displacement 
will vary approximately continuously with distance. We can regard it as a continuous function, 
x, with 

x(nAa) = x n 
The kinetic energy of the system is then simply 

1 N 
T=-AmY,x 2 n 
z o 

The total potential energy is similarly 

1 1 N 

V = -kx 2 Q + -kY^iXn- Xn-l) 2 
iL Z j 

From this we may deduce the equations of motion in two different ways. 
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6.4 Equations of motion: First approach 



If we take the large N limit first, these sums become integrals. 



lim Woo T = ^foX 2 ds 



where 5 is the distance from the wall, and 



= \kx 2 + \KZ ( ^ Aa+A A a )-^ Aa » f aAa 
lim Noo V = \kx{0) 2 + \Kaj{^ ds 

Since the springs always remain attached to the wall, x(0)=0. 

The integrand for T is the product of the density and the square of the velocity, just as we might 
naively expect. Simalarly, V is the integral of the potential energy of a infinitesimal spring over 
the length of the system. 

Using the Lagrangian will let us get equations of motion from these integrals. 
It is 



L = 



2a 



f a x 2 ds-\Kaf a (^fds 



The action of the system is 



fCC, m n Ka 2 „ 

/ Ldt= I I Se(x,x')dxdt wherei?= — [x z x' 2 ) 

J J J 2a m 

Here, we are integrating over both space and time, rather than just space, but this is still very 
similar in form to the action for a single particle. 

We can expect that the principle of least action will lead us to the natural extension of Lagrange's 
equations to this action, 

d d££ d gjf _ d<£ 
ds dx' dt dx dx 

This equation can be proven, using the calculus of variations. Using it for this particular La- 
grangian gives 

d 2 x Ka 2 d 2 x _ 
dt 2 m ds 2 
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This is a partial differential equation. We will not be discussing its solution in detail, but we will 
see that it describes waves. 

First, though, we will confirm that these equations are the same as we get if we find the equations 
of motion first, then take the large N limit. 

6.5 Equations of motion: Second approach 

First, we will look at the potential energy to see how it depends on the displacements. 

V = | IcXq + | k Y.q l {%n~ x n+i) 2 

= |fc(xQ ••• (x„_i - x n ) 2 + {x n - X n+ i) 2 ■■■ (X;v - Xjv-l) 2 ) 

= 2 (^Xq — 2xoXi 2x n x n+ \ + 2x^ — 2x n x n+ \ • • ■ x^ — 2xnXn-{) 

= 2 k {^ x ~ 2XoXi • • • 2x n [X n — X n + \ — X n -i) • • • X^ — 2XjvX]V_i) 

Notice that we need to treat the displacements of the first and last masses differently from the 
other coordinates, because the dependence of the Langragian on them is different. 

The kinetic energy is symmetric in the coordinates, 

1 N 
T=-AmY,x 2 n 
z o 

Using Lagrange's equations, we get that, for xq 

Amxo + 2fcxo = fcxi, 

for Xjv, 

Araxjv + fcxjv = fcxjv-i, 

and, for all the other x n 

Amx„ + 2kx n = k{x n+ i + x„_i), 
We can replace k and A m by the limiting values K and m using 

k Ka 2 
Am m\a 2 



Ka 2 Xi-2x 
m Aa 2 
Ka 2 x„_i+x„ + i-2x„ 

m Aa 2 
Ka 2 x N -i-x N 
m Aa 2 



giving us 

x 

X N 
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Looking at the general structure of the right hand side, we see difference between the displace- 
ment at nearby points divided by the distance between those points, so we expect that in the 
limit we will get differentials with respect to s, the distance from the wall. 

As N tends to infinity both xo and x\ tend to x(0), which is always zero, so the equation of 
motion for xo is always true. 

In the continuum limit, the equation for becomes 

Ka 2 1 dx 

x{a) = — — 

m Aa as 

Since A a tends to zero, for this to be true we must have x'=0 at a. 
For the other displacements, 

i: m x n -\+ x n+ i—2.x n i' 1 ( x n+ \—x n x n —x n - \ \ 

_ x'({n+\)Aa)-x'{nAa) 

~ llm AaO Aa 

_ d 2 x 



ds 2 



x-a 



so we get the equation 



d 2 x Ka 2 d 2 x _ 
dt 2 m ds 2 



just as with the other approach. 



6.6 Waves 



It is intuitively obvious that if we flick one of the masses in this system, vibrations will propagate 
down the springs, like waves, so we look for solutions of that form. 

A generic travelling wave is 



x = Asin(a)t -kx + a) 
Substituting this informed guess into the equation gives 

o K(l 2 n 

-Aa) sin(wf-KX + a) = k sin(aif- kx + a) 

m 

so this wave is a solution provide the frequency and wavenumber are related by 

o K& 2 •) 

oi z = k z 

m 

The speed of these waves, c, is 
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Thus, we've gone from Newton's laws to waves. 

We can do the same starting with a three-dimensional array of particles, and deduce the equa- 
tions for longitudinal and transverse waves in a solid. Everything we said about waves earlier will 
be true for these systems. 

This particular system has two boundary conditions: the displacement is zero at the wall, and a 
local extrema at the free end. This is typical of all such problems. 

When we take account of the boundary conditions we find that the correct solutions is a combi- 
nation of standing waves, of the form 



x = A m sin 



2b + Ins 



a 



<Ka 2 2b + 1 n 



sin 



t + a, 



m 



where b is any integer. 

If we also knew the initial displacement we could use Fourier series to obtain the exact solution 
for all time. 

In practice, N is typically large but finite, so the continuum limit is only approximately true. 
Allowing for this would give us a power series in 1/ N, describing small corrections to the ap- 
proximation. These corrections can produce interesting effects, including solutions, but we will 
not calculate them here. 



The continuum limit also fails for small wavelengths, comparable with the particle spacing. 



6.7 Fields 

In the continuum limit, the spring is described by a variable which is a function of both position 
and time. Variable such as this are commonly referred to as fields. 

At first sight, classical fields look quite different to classical particles. In one case position is the 
dependent variable; in the other, it is an independent variable. However, as the above calcula- 
tions suggest, fields and particles have an underlying unity, if we take a Lagrangian approach. 

We can deal with both using essentially the same mathematical techniques, extracting informa- 
tion about both the field and the particles in that field from the same source. 

E.g., once we know the Lagrangian for electromagnetism, we can deduce both the partial differ- 
ential equations for the EM fields, and the forces on charged particles in those fields, from it. We 
will see precisely how later, when we come to study electromagnetism. 

In the example above, the field Lagrangian was the continuum limit of a Lagrangian for the dis- 
crete system. It did not have to be. We can investigate the fields described by any Lagrangian we 
like, whether or not there is an underlying mechanical system. 
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So far, we've looked at waves and movement under Newton's law, and seen how the study of 
movement can lead us back to waves. Next, we will look at special relativity, and see how Ein- 
stein's insights affect all this. 
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